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We present a perturbative calculation of the renormalization constants for the AB = 2 four-quark operators, 
which is needed to evaluate the bag parameters Bb and Bs in lattice QCD. The one-loop coefficients are calculated 
for the 0(1/M) NRQCD heavy quark action and the 0(a)-improved Wilson action for light quark. Using these 
coefficients, which remove 0{as / aM) error, we reanalyze our previous calculation of the bag parameters. 



1. Introduction 

In the lattice calculation of the B meson B- 
parameters Q , the use of the effective theory for 
heavy quark is essential to obtain results with 
controlled systematic errors. The static approxi- 
mation in which the heavy quark mass is sent to 
infinity has been used by many authors and 
the next step is to incorporate the correction from 
finite h quark mass. In this direction, we per- 
formed an exploratory lattice calculation using 
the NRQCD action for heavy quark, in which we 
found that the l/ms correction could be sizable 
for Bb [§| and for Bs @]. A drawback in these 
calculations was, however, that the perturbative 
renormalization of the four-quark operators were 
not available for the NRQCD action and we used 
the one-loop coefficients for the static action. As 
a result our results contained large systematic er- 
ror of 0(as/aM), which is about 10% if evaluated 
with an order counting at (3 — 5.9. 

In this paper, we present a perturbative calcu- 
lations of renormalization constants of the heavy- 
light four quark operators with the 0(1/M) lat- 
tice NRQCD action. The results connect the 
missing link in our previous calculations, and we 
present a reanalysis oi Bb and Bs- A full de- 
scription of our results is published in [|| . 



2. Operator Definitions 

The B parameters in the Bq 
through 



Bq are defined 



{Bq\OrmB,) = ^MljlBf^ii,), 



(1) 



Bq\OT{^Ji)\Bq) 



Mf 



rubin) + mq{^) 



(2) 



where a subscript q denotes light quark flavor, 
i.e. c? or s. is usually called Bb, but we use 
a notation B^ in this paper to remind that it pa- 
rameterizes a matrix element of the operator Ol- 
In the above equations, the operators (/i) and 
0^^{^) are defined in the continuum theory with 
the MS scheme. Their Dirac structure is 



Ol = bjf_,PLqbji_,PLq, 
Os - bPLqlPLQ, 



(3) 
(4) 



*postor presented by K-I. Ishikawa. 



where Pl — 1 — 75. The totally anti-commuting 
convention for 75 is assumed in the dimensional 
regularization. 

The matching of these continuum operators is 
done against the lattice NRQCD. Namely we con- 
sider the 0{1/M) NRQCD action for heavy quark 
and the 0(a)-improved Wilson fermion for light 
quark. For gluons the standard plaquette action 
is used. 
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For Ol and Os the following heavy-light four- 
quark operators appear in the matching relation: 

Or = \hi.PRq\ [HPrq] , 

Os = \bPLq\ [bPLq] , 

On = 2 [b-f^PLq] [h^.PRq] + 4 [bPLq] [bPRq] , 
Om - 2 [h^.PLq] [h>.PB.q] - 4 [bPLq] [bPRq] , 
Op = 2 [h^.PLq] [h^.PRq] + 12 [bPLq] [bPRq] , 
Ot = {2 + TV,) [h,.PLq] [h^.PRq] 

-2{3N^ - 2N, - 4) [bPLq] [bPRq] , (5) 

where Pr = 1 + 75. The heavy quark field b 
has four Dirac components, and it is obtained 

from the two-component spinor non-rclativistic 
quark (anti-quark), Q (x), through the Foldy- 
Wouthuyscn-Tani transformation. 



1 



7 • D 
2M 



X' 



(6) 



In the static limit M 00, where alge- 
bra is greatly reduced, we consider the 0{a)- 
improvcment of the lattice four-quark operators 
through 0{aga). In doing that we also define the 
following dimension seven operators: 



Old 

OSD 



[b^^PLq] h^.PL{aD-^)q 



[bPLq] bPL{aD-j)q 



Ond = 2[b^^PLq 



O 



PD 



b^nPR{aD ■ ^)q 
+4 [6Pl9] ^PR{aD ■ 7)g 
2 [67^Pl9] [bl^,PR{aD ■ 7)9 
+12 [bPLq] \bPR{aD-^)q 



(7) 



3. Matching 



At one-loop level, the matching of the contin- 
uum operator 0^^(/i) is given by 



[ 



l + ^Cx,x(/x;l/a) 0'^\l/a) 
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Figure 1. Mass dependence of Cl,x- Ca is the 
matching coefficient for the temporal component 
of the heavy-light axial vector current. 



E 

Z=dim.7 



as 
An 



(:x,zOt{lla), 



(8) 



where the summations with Y and Z run over 
dimension six and seven operators respectively. 
fi is the rcnormalization scale in the MS scheme 
and a is lattice spacing. 

The coefficients Q are determined so that the 
on-shcU scattering amplitudes in the continuum 
and lattice theories agree with each other. For the 
external state, zero spatial momentum heavy and 
light quark suffice to determine the cocflicicnts for 
dimension six operators. A derivative in terms of 
spatial momenta have to be considered to obtain 
the dimension seven coefficients. 



r=dim.6 
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3.1. Ol 

The matching relation for Ol is 



CL,sil/a) Of{l/a) 



47r 



+^CL.M(i/a) o;r;(i/a) 

+ ^CL,«(l/a) Oi,^'(l/a) 



47r 



CL,LDOLD(l/a) 



(9) 



The coefficient (L.LifJ-', 1/a) of the leading opera- 
tor O^* contains logarithmic terms 61n(a^MQ) — 
21n(a^/x^), while others do not have fj, depen- 
dence. The operator Om does not exist in the 
matching of the static operator, but it becomes 
necessary for NRQCD. 

Figure 1 shows the heavy quark mass depen- 
dence of Cl,x for the dimension six operators. 
Filled symbols at l/aMg = are calculated for 
the static action. We find a large slope in 1/ uMq 
for as shown in the top panel of Figure 1. 
Most of the mass dependence comes from factor- 
ized diagrams, i.e. the four-quark operator can be 
split into two axial- vector currents and the gluon 
propagator does not connect these two. As a re- 
sult, the coefficient Cl,l is much reduced if we 
subtract the contributions of two axial vector cur- 
rent, as plotted in the lower panel (Cl,l — 2Ca). 
This combination will appear in the analysis of 
the B parameters, when we consider the ratio 

{OL)/{Aor. ^ 

The coefficients for the dimension seven oper- 
ators are obtained in the static limit as Cl,ld = 
-17.20 and Cl,nd = -9.20 §. 

3.2. Os 

The matching of Os is given by 

oTin) = [i + ^C5,5(/i;i/«)]or(i/a) 
+^C5,L(M;i/«)or(i/«) 

+ ^CsA^/a)0'^\l/a) 



40 



fi- s+41n(a2M6)-( 16/3)ln(a2^2) □ 
fi..5-2Ci+(i6/3)ln{M5//i2)+(8/3)ln(aW§) O 
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Figure 2. Mass dependence of Cs,x- 



+ ^CsM^/a)Ol^\l/a) 
+ ^Cs.SDOrDii/a) 
+ ^Cs.LDOTDil/a) 



+ ^Cs.PDO%il/a), 



(10) 



where Ot is new for the NRQCD action. The co- 
efficients Cs,x are plotted in Figure 2. Like Cl.l in 
the matching of Ol, the leading operator Os has 
large coefficient, which is reduced by subtracting 
the factorized contribution, as shown in the upper 
panel. 

The coefficients for the dimension seven opera- 
tors are Cs,SD = -6.88, Cs.ld = 2.58 and Cs,pd 
= 1.15 in the static limit. 
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4. Reanalysis of B Parameters 

With these perturbative coefficients, we rean- 
aiyze our previous iattice calculation for the B 
parameters, B^ and ,65 



4.1. Bl 

The B parameter B^ is calculated through 
B^{^^) =Y.^l.xiaA^: \la)B\\\la\ (11) 

X 

where X runs over L, S*, A^, M and R. The rcnor- 
malization constants are 



l + ^(CL.i(A*;l/a)-2a(l/a)) 



^L,SM2(l/a) 



477 



a,s(l/a), 



■^L,A'/A2(l/a) = ^C-L,Af(l/a)i 

■Z'L.M/A^Cl/a) = ^CL,M(l/a): 

^L,_R/A2(l/a) = ^CL,-R(l/a), 

and the lattice B parameters B^^^ (l/*^) were mea- 
sured in the numerical simulations through 



{B,\Of{l/a)\B,) 



l{B,\A't{^/o^mmT{^/o^)\B,y 

(12) 



in our previous paper ||3|,y. Since the dimen- 
sion seven operators are not included in the 
calculation, we could not remove the error of 

©(a^aAqcD)- 

Figure ^ shows the individual contributions 
jf/^2(/i; l/a)i?^*(l/a) for the B parameter 
Bl. The F-scheme coupling constant ay(g*) is 
used with q* = 2 /a. Although there is a siz- 
able mass dependence in x/A2(/^; 1/a) as we 
can see from the plot of coefficients (Figure 1), it 
is canceled by the 1/tob corrections in the ma- 
trix elements on the lattice i3l^*(l/a), and there 
is little mass dependence in the products shown 
in Figure ^. 

Our result for B^^ (fi) is plotted in Figure ^ as 
a function of 1/Mp with filled circles. The statis- 
tical error shown by solid error bars is very small 




1/MpiGeV ) 

Figure 3. Contribution of individual operators 
ZL.x/A^B'f{l/a) to BLiji), where X = L (O's), 
S (O's), R (>'s), N (A's), M (<'s). 



except for the heaviest heavy quark, while the sys- 
tematic error given by dotted error bars is more 
significant. We estimate the systematic errors 
with an order counting, in which the systematic 
errors of O(a^), O(asaAQCD), 0{asAQCB/M), 
0{{a Aqcb)^), and 0((Aqcd/M)2) are consid- 
ered and added in quadrature. For compari- 
son, we also show a result with one-loop coeffi- 
cients calculated for the static heavy quark by 
open circles, for which a large systematic error of 
0{as/iaM)) is added. The effect of the NRQCD 
renormalization constants is sizable, but still it is 
within the estimated systematic error. 

The numerical results for Bl at physical B me- 
son mass are 



Bl^m,) = 0.85(3)(11), 
Bf^imt) = 0.87(2)(11), 



(13) 
(14) 



where the first error is statistical and the second 
is systematic one 13%). The SU(3) breaking 
effect is 



^ = 1.01(1), 

Bf!{m,) 

where only the statistical error is quoted. 



(15) 
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Figure 4. 1/Mp dependence of i?||^^^(mb) (cir- 
cles) and B^^{m}j)/Tls{nibY (squares). Solid er- 
ror bar shows statistical error, dotted includes 
systematic. 

4.2. Bs 

Using Eq. (p^, Bs is obtained as 

B^U) ^ 

=Y,Zs,x/A^{t^; l/a)Bi^;(l/a), (16) 

where X runs over S*, L, R, P and T. In the LHS, 
we define a ratio of heavy-light matrix elements 



7^,(^) = 



(17) 



(0|PMS(^)|5^) 

The renormalization constants are 

Zs,s/A^{^J'■A/a) = 

l + ^(Cs,5(M;l/a)-2a(l/a)), 

^sxM2(Ai; 1/a) = ^Cs,L(l/a), 

^s,i?M2(l/a) = ^Cs,fl(l/a), 

Zs.p/A^C^/a) = ^Cs,p(l/a), 

"-Cs,T(l/a). 



47r ■ 



^s,TM2(l/a) 
For we define the lattice B parameters 



as 



{B,\Ofil/a)\B,) 
-l{B,\Ar{l/a)\0){0\Aril/a)\B,)- ^ ^ 



The mass dependence of Bg^^/TZg is plotted 
in Figure ^ (filled squares). Open squares are 
our previous results calculated with the static 
renormalization constants. The mass dependence 
is less significant in the results with the correct 
(NRQCD) renormalization constants. 

The numerical results are 



J) d = 1.24(3)(I6), 



Bf^irrib) 



and 

BYf{mb)na{mbf 
Bff{mb)ns{mbY 



1.24(3)(I6), 



1.003(4). 



(19) 
(20) 

(21) 



for the SU(3) breaking ratio. 
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